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I. Summary from previous reporting periods 


This is the final progress report covering the complete period from December 1, 1988 to 
August 31, 1991, funded under the NASA Contract NAG- 1-922. The research effort was 
directed toward the problems of real-time trajectory optimization and guidance law 
development for National Aerospace Plane applications. In particular, singular perturbation 
methods were used to develop guidance algorithms suitable for on board, real-time 
implementation. 

The bulk of the work completed during the period is summarized in three previous 
reports, listed as references [l]-[3]. Reference [4] cites the final progress report from the 
previous project (funded under the NASA Contract NAG- 1-784). 

The work completed during the period from December 1, 1988 to June 30, 1989 (See 
Ref. [1]) consisted primarily of extensions to the analysis reported in reference [4] to 
include a number of important considerations. In particular, the vehicle model was 
extended to include angle of attack effects, the thrust vector component normal to the 
velocity vector, and flight in the subsonic and supersonic regimes. A multi-mode 
propulsion system consisting of turbojet, ramjet, scramjet and rocket engines was assumed 
and simple models for thrust generation and fuel consumption were adopted for each 
engine cycle. The state-space was further constrained by considering a maximum allowable 
heating rate. Singular perturbation methods were applied to this more realistic model, 
leading to a simple algorithm suitable for generating a nearly-fuel-optimal altitude profile in 
real time. A simple iterative algorithm was derived that approximates the optimal engine 
transition points and the regions of cycle overlap. Feedback linearization was employed to 
derive an angle of attack controller which can be used to guide the vehicle along the nearly- 
fuel-optimal altitude profile in simulations of flight within the atmosphere. A computer 
subroutine based on the space shuttle explicit guidance algorithm was written to handle the 
exoatmospheric phase of ascent guidance which allows for the simulation of insertion into 
orbit. The resulting software was employed to examine the influence of the added model 
complexity on the fuel-optimal ascent trajectories and the performance of the guidance 
algorithms^ 

During the second reporting period, up to December 31, 1989 (See Ref. [2]) general 
problems associated with on-board trajectory optimization, propulsion system cycle 
selection, and with the synthesis of guidance laws were addressed for ascent to low-Earth- 
orbit of an air-breathing, single-stage-to-orbit vehicle. The work built directly upon the 
analytical results of reference [1]. A good portion of the work focused on making 
improvements to the vehicle models employed. The NASA "Generic Hypersonic 
Aerodynamic Model Example" and the "Langley Accelerator" aerodynamic data sets were 
acquired and implemented. Work pertaining to the development of purely analytic 
aerodynamic models also continued at a low level. A generic model of a multi-mode 
propulsion system was developed that includes turbojet, ramjet, scramjet, and rocket engine 
cycles. Provisions were made in the dynamic model for a component of thrust normal to 
the flight path. Computational results, which characterize the nonlinear sensitivity of 
scramjet performance to changes in vehicle angle of attack, were obtained and incorporated 
into the engine model. Additional trajectory constraints were also introduced. The 
constraints treated were: maximum dynamic pressure, maximum aerodynamic heating rate 
per unit area, angle of attack and lift limits, and limits on acceleration both along and 
normal to the flight path. 

The remainder of the research effort during the second period focused, for the most part, 
on required modifications to the previously derived algorithm when the model complexity 
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cited above was added. In particular, analytic switching conditions were derived which, 
under appropriate assumptions, govern optimal transition from one propulsion mode to 
another for two cases: the case in which engine cycle operations can overlap, and the case 
in which engine cycle operations are mutually exclusive. The resulting guidance algorithm 
was implemented in software and exercised extensively. It was found that the 
approximations associated with the assumed time scale separation employed in this work 
are reasonable except over the Mach range from roughly 5 to 8. This phenomenon is due to 
the very large thrust capability of scramjets in this Mach regime when sized to meet the 
requirement for ascent to orbit. Very little mass penalty is induced by the resulting 
inaccuracies in the trajectory over this region because it is traversed rapidly. However, the 
reduced solution climb paths prove to be unfeasible within this Mach range when subject to 
the full model dynamics and active trajectory constraints. These difficulties were 
successfully overcome by accounting for flight path angle and flight path angle rate in 
construction of the flight path over this Mach range. The resulting algorithm provides the 
means for rapid near-optimal trajectory generation and propulsion cycle selection over the 
entire Mach range from take-off to orbit given a realistic nonlinear vehicle model and all 
pertinent trajectory constraints. 

A significant problem area encountered was the lack of a general theory for singularly 
perturbed systems that are subject to state-variable inequality constraints (Ref. [2]). Such 
constraints are common to a wide class of flight vehicles but have received little attention in 
the literature when the dynamic system is singularly perturbed. A study was initiated in this 
area and it was found that, when the reduced solution lies on a state-variable inequality 
constraint boundary, the boundary layer trajectories are of finite time in the stretched time 
scale. The possibility of costate discontinuities at the juncture between constrained and 
unconstrained arcs makes direct application of existing theory difficult at best. A 
transformation technique was identified that eliminates some of these difficulties, but at the 
cost of possibly increased system order and the introduction of singular arcs. Much work 
remains to be done in this area. 

Work on the development of simple, efficient algorithms for prediction of vehicle 
aerodynamic and propulsive performance continued during the second phase of the 
program (Ref. [2]). Improvements in the modeling of the hypersonic lifting body module 
eliminated previous discrepancies between measured and predicted aerodynamic behavior. 
Several modes of data entry can now be implemented making assessment of a given vehicle 
configuration very simple. An interactive program mode was devised that makes possible 
direct and immediate assessment of configuration changes on selected vehicle performance 
parameters. The algorithms developed in this program are of potential use in applications 
beyond those originally envisioned. 

The first two reporting periods resulted in four conference papers (Ref. [5]-[8]) which 
discuss most of the results of this research effort. A Ph.D. Dissertation that details the 
entire effort to date was published in December of 1989 (Ref. [9]). A full-length paper 
entitled "Rapid Near-Optimal Trajectory Generation for Single-Stage-to-Orbit Airbreathing 
Vehicles" has also been submitted and accepted for publication in the AIAA Journal of 
Guidance, Control and Dynamics (Ref. [10]). 

There was a funding lag for the period from January to May of 1990 during which no 
research was conducted. Funding resumed on the first of May. 

During the next reporting period, from May 1 to October 31, 1990 efforts were primarily 
focused upon developing a general understanding of singularly perturbed systems subject 
to state-variable inequality constraints and also upon developing criteria for the applicability 
of singular perturbation techniques to flight mechanics problems and in particular to aircraft 
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energy climbs (Ref. [3]). As noted in reference [2], singularly perturbed optimal control 
problems with state-variable inequality constraints can exhibit complex boundary layer 
phenomena. In particular, the boundary layer transitions associated with such problems 
can be of finite time when the state constraint is first encountered at the end of the boundary 
layer transition. The lack of a general theory for treating such systems was identified as a 
significant research problem. 

A cursory look at the problem was completed prior to the submission of the 1989 final 
report (Ref. [2]). Since that time, considerable progress was made. The results of the 
effort are detailed in a technical paper that was first presented at the 1990 AIAA GN&C 
Conference (Ref. [1 1]). A revision of the paper has also been accepted for publication in 
the AIAA Journal of Guidance, Dynamics and Control. The results are summarized as 
follows. 

The established necessary conditions for optimality in nonlinear control problems that 
involve state-variable inequality constraints were applied to a class of singularly perturbed 
systems. The distinguishing feature of this class of two-time-scale systems is a 
transformation of the state-variable inequality constraint, present in the full order problem, 
to a constraint involving states and controls in the reduced problem. The existence of a 
nonsingular control solution was assumed. It was of particular interest to construct the 
zeroth order initial boundary layer solution, or at least and approximation to it, when the 
reduced solution lies on a state constrain boundary. In the absence of a state constraint, 
one can take advantage of the fact that the reduced solution serves as an equilibrium point 
for the boundary layer system. However, it was shown that, when a state constraint is 
active in the reduced problem, the boundary layer problem can be of finite time in the 
stretched time variable. Thus, the usual requirement for asymptotic stability of the 
boundary layer system is not applicable, and can not be used to construct approximate 
boundary layer solutions. Furthermore, an active state constraint introduced die possibility 
of discontinuous costate variables at the juncture between constrained and unconstrained 
arcs. 

Various means for treating such problems were investigated. A simple linear example 
was constructed and used to show that a Valentine transformation can be used to regain 
smoothness, but with limited advantage. That is, Valentine's transformation can be used to 
avoid the problems associated with discontinuous costate time histories, but at the expense 
of introducing a singular arc and discontinuities in the transformed control variable. A 
second linear example was used to illustrate the exact analytic solution of a simple 
singularly perturbed problem involving a state variable inequality constraint. The solution 
includes a "fast" costate discontinuity and a finite-time initial boundary layer transition. A 
third, but nonlinear, example for which the boundary layer system could not be solved 
analytically was then constructed. This example was used to illustrate a general feedback 
strategy that was developed for synthesizing a near-optimal boundary layer transition onto a 
constrained arc. In this technique, the costate jumps that can occur and the boundary layer 
final time are used as free parameters in order to satisfy continuity conditions in the state 
variables at the end of the boundary layer response. The resulting approximation was 
compared directly with the numerically generated optimal solution. The method proved 
quite satisfactory when used to construct an approximate solution for this relatively simple 
nonlinear example, at least for small perturbations away from the reduced solution. 

Several problems requiring further attention were identified. For instance, numerical 
problems were sometimes encountered in the solution process for Example 3 (See Ref. 
[11]) as time-to-go approached zero (i.e. as the boundary layer transition nears 
completion). This difficulty did not prevent the generation of an accurate approximation of 
the optimal solution for the example problem, and further manipulation of Example 3 has 
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lead to a completely analytic characterization of the solution. However, the possibility of 
approaching a singularity should be investigated in a generic setting. The approximation 
technique being employed depends upon a linearization of the boundary layer system about 
a non equilibrium point. It does not appear possible to characterize the stability of the 
approximation for a given set of initial conditions. And though guidance along a constraint 
boundary will likely be subject to small perturbations only, linearization does introduce the 
likelihood of control saturation for sufficiently large perturbations. A multiple time scale 
approach, in which altitude and flight path angle dynamics are examined on separate time 
scales, could eliminate this dependence on linearization. There is also a question of 
applicability when atmospheric disturbances lead to a constraint violation. 


II. P rogress this reporting period 

During the last reporting period, from November 1, 1990 to August 31, 1991 efforts 
were again focused upon developing a general understanding of singularly perturbed 
systems subject to state-variable inequality constraints and also upon developing more 
stringent criteria for the applicability of singular perturbation techniques to flight mechanics 
problems and in particular to aircraft energy climbs. Specifically, a systematic approach 
was devised for naturally identifying the perturbation parameter e in a singular perturbation 
analysis of aircraft optimal guidance and expressing it in terms of original physical problem 
parameters. The approach, which is based on a nondimensionalization of the equations of 
motion, can be used to evaluate the appropriateness of forced singular perturbation 
formulations used in the past for transport and fighter aircraft, and to assess the 
applicability of energy state approximations and singular perturbation analyses for 
airbreathing transatmospheric vehicles with hypersonic cruise and orbital capabilities. 
Furthermore, the approach can easily be extended to asses the possibility of treating the 
same problems by assuming multiple (more than two) time scale behavior. 

The lack of strict criteria for the applicability of singular perturbation techniques to flight 
mechanics problems served as an incentive for our efforts during the last reporting period. 
The methods of matched asymptotic analysis in singular perturbation theory are based on 
the presence of small parameters in the differential equations of motion which give rise to 
multiple time scale behavior. It has been noted by numerous authors (Ref. [13], [14]) that, 
in spite of a wide number of papers attesting to the applicability of singular perturbation 
methods to optimization problems in aircraft flight mechanics, few have been successful in 
first casting the equations of motion in a singular perturbation form. A few notable 
exceptions are Refs. [13]-[16]. In Ref. [13] two methods for time scale separation analysis 
are proposed to identify the proper assignment of state variables to various time scales. 
These methods are based on forming an estimate of the state variable speeds. In Ref. [14] a 
rescaling to nondimensional variables is recommended. However, it is noted that the proper 
scaling transformation is not obvious, even if the time scale separation of the variables is 
well known from analysis or experience. Both of these papers (and in particular Ref. [13]) 
provide extensive references to earlier studies which employ so-called forced singular 
perturbation formulations, in which the perturbation parameter (say e), nominally equal to 
1, is artificially introduced as a book keeping parameter in a formal expansion of the 
solution about e = 0. In particular, there exists a large number of publications on the 
optimization of aircraft energy climbs (see for example Refs. [17]-[20]), none of which 
make any attempt to identify an appropriate perturbation parameter in terms of the original 
problem parameters. This is particularly disturbing considering the number of years that 
have passed since such analysis techniques were first introduced in the flight mechanics 
literature. In any singular perturbation analysis, every attempt should be made to identify 
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the perturbation parameter in terms of the original problem parameters (which in general 
include the boundary conditions) so that the physical process that gives rise to the two time 
scale behavior is clearly understood. Then, the range of parameter values for which the 
perturbation analysis is valid can be easily identified. Knowledge of time settle separability 
present in the system dynamics, and success in exploiting this characteristic to obtain 
approximate solutions, is not, in itself justification for artificially introducing e. That is, 
within the framework of our system of logic it is always possible to have conclusions that 
are true , which follow from assumptions that are wrong. 

Our work has partially rectified this situation by presenting a systematic (albeit still ad- 
hoc) approach to nondimensionalize variables in nonlinear optimization problems in flight 
mechanics. Most of the considerations that were presented apply in other fields as well. 
Our main motivation for collecting and stating these considerations was to define the 
thought process by which it is possible to arrive at a suitable scaling of the aircraft energy 
climb problem. Of particular interest was the assessment of the applicability of energy state 
approximations and singular perturbation analyses for airbreathing transatmospheric 
vehicles with hypersonic cruise and orbital capabilities. 

The major result of our effort was the demonstration that for energy climbs that take 
place on a vertical plane the singular perturbation parameter e is always equal to the 
maximum longitudinal loading factor of the vehicle. Two time scale behavior is suggested 
according to whether e is less than one or not. Based on this result it is straightforward to 
see why singular perturbation methods applied to aircraft performance optimization have 
worked so well in the past. The maximum longitudinal loading factors associated with the 
majority of conventional aircraft are either less than one because the aircraft lack very high 
thrusting capabilities, or because they are restricted to be so for other reasons (structural, 
comforting, etc.). A few notable exceptions do occur for some modem fighters. This 
directly suggests that most conventional aircraft can be expected to exhibit two-time- scale 
behavior for almost any energy climb that they are allowed to perform. This then appears to 
be the reason for the past success of so many singular perturbation treatments of aircraft 
energy climbs. The implication for transatmospheric vehicles is rather direct. If we consider 
such a vehicle as a passenger transport, then, in order to assure passenger comfort it is only 
natural to impose as a constraint a maximum longitudinal loading factor for the vehicle that 
is less than one. Our work suggests that such a constraint would imply two-time-scale 
behavior for any type of energy climb that such a vehicle would be allowed to perform. 
Therefore, singular perturbation formulations of such maneuvers still appear to be 
promising, even with all the added complexities that the flight regimes of such vehicles can 
involve. A paper which reports on our progress has been presented at the 1991 AIAA 
GN&C Conference. A copy of the paper is included as an Appendix. 
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Abstract 


This paper proposes a systematic approach for identifying the perturbation parameter in singular 
perturbation analysis of aircraft optimal guidance, and in particular considers a family of problems 
related to aircraft energy climbs. It is first shown that for energy climbs that take place on a vertical plane 
the singular perturbation parameter can always be taken to be the maximum allowed longitudinal loading 
factor of the vehicle. Two time scale behavior is suggested according to whether this parameter is 
sufficiently less than one. The approach, which is based on a nondimensionalization of the equations of 
motion, is then used to evaluate the appropriateness of forced singular perturbation formulations used in 
the past for transport and fighter aircraft, and to assess the applicability of energy state approximations 
and singular perturbation analysis for airbreathing transatmospheric vehicles with hypersonic cruise and 
orbital capabilities. 


Introduction 

The methods of matched asymptotic analysis in singular perturbation theory are based on the 
presence of small parameters in the differential equations of motion which give rise to multiple time scale 
behavior. It has been noted by several authors 1 - 2 that, in spite of a wide number of papers attesting to the 
applicability of singular perturbation methods to optimization problems in aircraft flight mechanics, few 
have been successful in first casting the equations of motion in a singular perturbation form. A few 
notable exceptions are Refs. 1-4. In Ref. 1 two methods for time scale separation analysis are proposed 
to identify the proper assignment of state variables to various time scales. These methods are based on 
forming an estimate of the state variable speeds. In Ref. 2 a rescaling to nondimensional variables is 
recommended. However, it is noted that the proper scaling transformation is not obvious, even if the 
time scale separation of the variables is well known from analysis or experience. Both of these papers 
(and in particular Ref. 1) provide extensive references to earlier studies which employ so-called forced 
singular perturbation formulations, in which the perturbation parameter (say e), nominally equal to 1, is 
artificially introduced as a book keeping parameter in a formal expansion of the solution about £ = 0. In 
particular, there exists a large number of publications on the optimization of aircraft energy climbs (see 
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for example References 5-8), none of which make any attempt to identify an appropriate perturbation 
parameter in terms of the original problem parameters. This is particularly disturbing considering the 
number of years that have passed since such analysis techniques were first introduced in the flight 
mechanics literature. In any singular perturbation analysis, every attempt should be made to identify the 
perturbation parameter in terms of the original problem parameters (which in general include the 
boundary conditions) so that the physical process that gives rise to the two time scale behavior is clearly 
understood. Then, the range of parameter values for which the perturbation analysis is valid can be 
easily identified. Knowledge of time scale separability present in the system dynamics, and success in 
exploiting this characteristic to obtain approximate solutions, is not, in itself justification for artificially 
introducing e. That is, within the framework of our system of logic it is always possible to have 
conclusions that are true, which follow from assumptions that are wrong. 

In this note we attempt to partially rectify this situation by presenting a systematic (albeit still ad- 
hoc) approach to nondimensionalize variables in nonlinear optimization problems in flight mechanics. 
Most of the considerations that are presented apply in other fields as well. Our main motivation for 
collecting and stating these considerations is to define the thought process by which it is possible to 
arrive at a suitable scaling of the aircraft energy climb problem. Of particular interest is an assessment of 
the applicability of energy state approximations and singular perturbation analysis for airbreathing 
transatmospheric vehicles with hypersonic cruise and orbital capabilities. 


Subsonic-Supersonic Regimes, Flat Earth Approximation 
Consider atmospheric flight of a conventional aircraft, viewed as a point mass, in a vertical plane 
over a flat Earth. The equations governing such flight can be reduced to a three-state model in: mass 
specific energy E, flight path angle y, and altitude h. The vehicle mass, m, is assumed to be constant. 
The equations are: 


dE _ V(T - D) 
dt m 

dy ( L N ( g cos y\ 
dt ~ v mV ) \ V ) 


( 1 ) 


1 1 


( 2 ) 


= V sin y 


(3) 


where L, D and g denote the lift, the drag and the (constant) gravitational acceleration. It is assumed that 
the atmosphere is stationary, and that the thrust, T, is directed along the flight path. The specific energy 
(mechanical energy per unit mass of the vehicle) E and the speed V are related by 



(4) 


and E rather than V has been employed as a state variable. 

In many of the earlier singular perturbation studies the traditional way of writing down Eqs. (2) and 
(3) was: 

dy ( l (Z cos y\ 

e dT = lmVV l V J ( 5 ) 

= v sin y (6) 

that is, by artificially introducing a parameter e and then stating that its nominal value was equal to 1. 
Since our main purpose in the present paper is to avoid such an artificial introduction at the outset, Eqs. 
(5) and (6) will serve only as a guide for the natural introduction of e. 


Nondimensional form 

The first step is to put Eqs. (1) through (3) in nondimensional form. To this end we define the set S 

S = {t 0 , E 0 , h Q , Vq, T q , Dq, ( 7 ) 

The elements of the set S are at this point arbitrary positive quantities, and the only restriction that we 
impose upon them is that: 

to has dimensions of time 
E 0 has dimensions of energy per unit mass 
ho has dimensions of lenght 
V 0 has dimensions of speed 
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T 0 , D 0 , and Lo have dimensions of force 
Using the elements of S to define the nondimensional quantities: 



Eqs. (1) through (3) can be put into the following nondimensional form: 


_ vr it - d d ) { t() V ° 1 
dt V ('T 0 UL V^E 0 m J 

(10) 

dy i 

dt” - ' 

^LY L o t o > l rcosYYg^"! 
Lv^mV 0 J A V A V oJ 

(ID 

dh 

dt 

fW t \ 

V sin Y 

h o J 

(12) 


The goal is now to put Eqs. (10) through (12) in the traditional singular perturbation format. We thus 
multiply both sides of Eqs. (1 1) and (12) by (h 0 / V 0 to ). This results in: 

' h o VY _ f L Y L o h o l _ f cos V \( g V 

IVoM IvJ mV 2 l V J v 2 (13) 


. V 0 l 0 J 


dh 

dt 


= V sin y 


(14) 


Comparing the set of Eqs. (10), (13) and (14) with the set (1), (5) and (6), it is evident that we can make 
the two sets similar by imposing the following four conditions on the elements of the set S: 

(15) 

(16) 

(17). 


T = D 
A o o 

T t V 
1 0 1 0 v 0 

E 0 m 


= 1 


L o h o 

mV„ 2 


= 1 
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( 18 ) 



= 1 


If we define e as 


then, Eqs. (10), (13) and (14) assume the form: 

= V(T-D) 

dy ( L - cos y'] 

e dT = l V J 

dh . 
e— = V sm y 
dt 


(19) 


( 20 ) 

( 21 ) 

( 22 ) 


To summarize, it was shown in the present section that it is possible to introduce a parameter e 
naturally into the equations of motion (Eqs. (1) - (3)) by first introducing a set of arbitrary positive 
quantities S (see Eq. (7)) to scale the variables of interest, and then by imposing four conditions (Eqs. 
(15) - (18)) on these quantities so that the resulting nondimensional equations assume the traditional 
singular perturbation form (Eqs. (20) - (22)). Note that only one of the arbitrary quantities in S is 
uniquely determined at this point. Combining Eqs. (17) and (18) it follows that 

L o = m g (23) 


Specifying a particular nondimensional form 

As shown in the previous section, only /our conditions are imposed on the seven elements of the 
set S in transforming the equations of motion to the traditional singular perturbation format. This means 
that we can specify three of the elements of S to fit our convenience and then determine the remaining 
four using Eqs. (15)-(18). The first conclusion therefore is that in general the value of e is quite 
arbitrary. For example, by choosing h 0 , Vo and to in two different ways e can be made arbitrarily small 
or large. The separability of the time scales on the other hand is a property of the system and not of the 
particular nondimensional form of the equations of motion that is chosen. We therefore expect that if the 
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system does indeed have this property it will exhibit it no matter what the actual value of e is. This is 

precisely the reason for the success of so many singular perturbation treatments of the past in which e 

was introduced artificially and its nominal value was said to be fixed at one. 

Although there is no unique way of specifying a particular nondimensional form of the equations of 

motion, we will now argue that there is one that results in additional physical insight. First, in order to 

maintain the relationship in Eq. (4) in the transformed variables, a fifth condition is introduced 

E„=gh„ (24) 


which together with Eq. (18) gives 



Using Eqs. (16) and (24) in Eq. (19) it follows that 

e= A 

fc mg 


(25) 


(26) 


It is now evident that only two among the seven elements of the set S need to be specified. Then, the 
five conditions, Eqs. (15) - (18) and (24), uniquely determine the remaining elements. 

Eq. (26) implies that e depends only on T 0 and is independent of the value of the remaining 
elements of S. The question therefore arises naturally as to whether there is a particular choice of T 0 for 
which the resulting value of e can be used as a strict criterion for the applicability of a singular 
perturbation analysis to Eqs. (20)-(22). The answer to this question is negative because, in a given time 

interval, it is the relative magnitudes of the three quantities 

dE dy dh 
dt ’ dt ’ dt 


and the boundary conditions of interest that determine the validity of a singular perturbation analysis. 
Specifically, for an aircraft to exhibit the well-known two-time-scale behavior in a given time interval, it 
is necessary that: 


dE dy 
dt ' ' dt 


(27) 
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dE / . dh 
dt U dt 


( 28 ) 


in that interval, and that the required change in E is sufficiently large to permit the boundary layer 
responses in h and y to reach their equilibrium values. Hence, the very question is whether or not the net 
change in E during the boundary layer response is sufficiently small to permit approximating E as a 
constant (to zero order in e) in the boundary layer analysis. In addition, we are interested in knowing if 
this two-time- scale property is a consequence of the inherent dynamics of the aircraft, and not a 
consequence of using a high gain control solution for L Therefore, we assume that the L resulting from 
the boundary layer analysis is of order one in Eq. (21). 

Under the above assumptions, there is a choice for T 0 for which the value of e can be used as a 
measure for the existence of time intervals in which two-time-scale behavior is exhibited. If the choice 
of T 0 is such that dE/dt is at most of the same order of magnitude as edy/dt and edh/dt, then, a value of 
e sufficiently less than one indicates the possible existence of such intervals. By suitably choosing V 0 we 

can restrict V to be of order one. Then, for the choice 

T 0 = ( T - D U (29) 


dE/dt is of order one, and both dy/dt and dh/dt are of order 1/e. For this choice of T 0 , e is given by 

(T-D) 

' max 

e = mg (30) 


and is equal to the maximum longitudinal loading factor of the vehicle. 

Note that Eq. (30) actually represents an upper bound for e since it is obtained by selecting the 
flight condition where the difference between thrust and drag, T-D, reaches a maximum. The logical 
choice for V 0 is the speed at this flight condition. One can also adopt the viewpoint of evaluating e along 
the energy climb path that results from a reduced solution. The value of e as a function of E can then be 
used as a measure to distinguish energy levels where a singular perturbation analysis may be appropriate 
from other levels where it may not be valid. 

It is interesting to note that a good deal can be anticipated from Eq. (30) for conventional aircraft 
without exact numerical evaluation: (T-D) max divided by mg is approximately equal to siny max where 
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y max is the maximum climb angle that can be maintained at a given energy level without loss of airspeed. 
It follows therefore that eel for all such aircraft types. For transport aircraft siny max is approximately 
0.1, while for fighter aircraft siny max is approximately 0.8. This suggests that the forced singular 
perturbation analysis used in the past studies of optimal aircraft trajectories is valid for most conventional 
subsonic and supersonic aircraft. 

Eq. (30) can also be used to estimate e in terms of the quantities (T/mg),,^ and (L/D) max for a given 
aircraft. Since L is less than or equal to mg along the energy climb path, it follows that a second upper 
bound for e is given by 

e<£uB (31) 

where eyg is defined as 

e UB s [( T/m S> m „~ 1/(L/D) m „] (32) 


Estimates of are given in Table 1. 

Table 1 

Estimation of based on Eq. (32) 


Parameter Transports Fighters 

TTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTvTTT 1 r*i* , *r 


(T/mg) 0.25 0.90 

max 

(WD) max 13-15 4-7 

e UB 0.17 - 0.18 0.65 - 0.76 


Hypersonic Regime 

Consider the flight of a hypersonic and possibly transatmospheric vehicle, viewed as a point mass, 
in a vertical plane over a spherical non-rotating Earth. The equations governing such flight can be 
reduced to a four-state model in E, m, y and radial distance from the center of the Earth, r. The equations 


arc: 

dE V(qT-D) 
dt ” m 


(33) 


1 7 




( 34 ) 


dy 
d T 



dr 

dt 


= V sin 


Y 


(36) 


where T is the maximum available thrust at a given speed and altitude. The control variables are L and t), 


where 0<rj<l is introduced as a nondimensional throttling variable. E and V are now related by 

,.2 


E = 


£ 

r 


(37) 


Note again that in earlier singular perturbation studies^, the traditional way of writing down Eqs. 
(35) and (36) was: 




= V sin 


Y 


( 38 ) 

( 39 ) 


that is, by artificially introducing a parameter e and then stating that its nominal value was equal to 1. 

Again, in order to avoid this artificial introduction, Eqs. (38) and (39) will serve only as a guide for the 
natural introduction of £. 


Nondimensional Form 

In order to put Eqs. (33) through (36) in nondimensional form we now define the set of arbitrary 
positive quantities 

S-{to.Eo,m 0 ,ro,V 0 ,fo,T 0 ,D 0 ,Lo} (40) 

and impose the restrictions that: 
to has dimensions of time 
E 0 has dimensions of energy per unit mass 
mo has dimensions of mass 
r 0 has dimensions of lenght 
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V 0 has dimensions of speed 
f 0 has dimensions of mass per unit time 
T 0 , D 0 , and Lo have dimensions of force 
Using the elements of S to define the nondimensional quantities: 


t = — — • E= — • m = - m - • r= ’ V = 
1 t„ » c E. ' 1,1 m n ’ r r n » v 


f=X. j=— • D=— • L= — 

f o T o D o ’ L o 

Eqs. (33) through (36) can be put into the following nondimensional form: 
dE V(ilTT 0 - DD 0 )f t 0 V 0 ^ 


(41) 

(42) 


dt 

dm 

dt 


m 


V E o ra o. 


(Ho 

l m o 


f(r, V, tj) 


dy 

f-M 

Wo ) 

/'cos y\ 

f »“o Y fVc OS yVVo^") 

d T = 

VmV J 

l m o v oj 

' l Vr 2 J 

[v 0 r 2 J A r A r o J 




(43) 

(44) 

(45) 

(46) 


In order to put Eqs. (43) through (46) in the traditional singular perturbation format, we multiply both 
sides of Eqs. (45) and (46) by (r 0 / V 0 to ). This results in: 

r o Yy f l V L o r o V f cosy y n 
Am v ^m 0 V 2 0 J'lvr 2 J 

( r ° Yr = y 


y 2 r 
V v 0 r 0. 




cos y 


(47) 

(48) 


IVoJdt 

Comparing the set of Eqs. (43), (44), (47) and (48) with the set (33), (34), (38) and (39) results in the 
following five conditions on the elements of the set S: 

T o =D o (49) 
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( 50 ) 


T t V 

A o l o v o 

E o m o 


= 1 


m „ 


= 1 


(51) 


L o r o 

■”o V o 


= 1 


(52) 




v o r o 


= 1 


(53) 


By defining 


e = 


v o^o 


(54) 


Eqs. (43), (44), (47) and (48) assume the traditional singular perturbation form: 

dE _ y(n T~ D) 
dt 


m 


dm 

dt 


= - f(r, V, ri) 


(55) 

(56) 


dy 

f L) 

/'cos y\ 

(V cos y^ 


e dF 

Imvi * 1 

l Vr 2 r 

l r J 

(57) 

dt 

= V sin y 



(58) 


Specifying a particular nondimensional form 

For the hypersonic case only five conditions on the nine elements of the set S are needed in order 
to put the equations of motion in the traditional singular perturbation format Thus, we can specify four 
of the elements of S to fit our convenience and then determine the remaining five using Eqs. (49) - (53). 

Again, in order to maintain the relationship in Eq. (37) in the transformed variables, a sixth 
condition is introduced 
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(59) 


E 


o ~ 


iL 

r o 


which together with Eq. (53) gives 



(60) 


If we think of r 0 as a radial distance, then Eq. (53) restricts V 0 to be the circular orbital speed at r 0 . 
Similarly, Eq. (52) restricts L 0 to be the centrifugal force that a point mass mo would experience in a 


circular orbit at r 0 . Using Eqs. (50), (53), and (59) we have 

T r 2 
o 

e “ 


(61) 


Hence, by picking three among the nine elements of the set S arbitrarily, the six conditions Eqs. (49) - 
(53), and (59) uniquely determine the remaining elements. 

The question arises again as to whether there is a particular choice for these three elements for which 
the resulting value of £ can be used as a measure for the applicability of a singular perturbation analysis 
to Eqs. (55)-(58). The right-hand- side of Eqs. (55) and (58) can be made of the same order of 
magnitude by choosing the ratio T 0 /mo as 



(62) 


Choosing r 0 as the sea level radius r SL , r and V are of order one. Also, for these choices of T 0 /mg and 
r 0 , dE/dt is of order one, and both dy/dt and dr/dt are of order 1/e. By choosing f 0 as the value of f at 
the flight condition where the ratio (T]T- D)/m is a maximum, dm/dt can also be made of order one. 


With the above choices of T 0 /mg, r 0 , and f 0 

r 2 

e = ^ 

l M- J 


(^) 


(63) 


max 


The right-hand-side of Eq. (63) is the maximum longitudinal loading factor of the vehicle in units of 
sea-level g's, and actually represents an upper bound for e since it is obtained by selecting the flight 
condition where (rjT- D)/m reaches a maximum. One can again adopt the viewpoint of evaluating e along 
the energy climb path that results from the reduced solution. The value of £ as a function of E can then be 
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used as a measure to distinguish energy levels where a singular perturbation analysis may be 
appropriate, from other levels where it may not be valid. 

A hypersonic flight vehicle employing an airbreathing propulsion system and sized for acceleration 
to orbital velocity necessarily employs a multimode propulsion system. An example might include 
turbojet, ramjet, scramjet, and rocket modes. Each mode of propulsion can be characterized by a 
corresponding e. Current models of this vehicle type exhibit large values of excess thrust at low 
hypersonic Mach numbers. In fact, Eq. (63) will produce an e that is greater than one over such flight 
phases. Experience with hypersonic vehicle dynamics reported in Ref. 10 indeed suggests that the 
assumed time scale separation is not valid in these phases. However, over the majority of the trajectory, 
Eq. (63) results in an e that is less than one just as in the Flat Earth, Subsonic-Supersonic case. 

Numerical validation 

It was shown in the preceding sections that for aircraft energy climbs that take place in a vertical 
plane, the singular perturbation parameter e can always be identified as the maximum longitudinal 
loading factor of the vehicle, measured in units of sea-level g's. In order to further explore the 
implications of this result, numerical evaluations of e will be presented in this section for several types of 

vehicles. 

The idea that the authors would like to introduce at this point is that in general, for a given aircraft, it 
may be sufficient to evaluate an upper bound for e, valid for the entire envelope, in order to get a hint 
for the possible two-time-scale behavior of the aircraft in question. If the resulting value of this upper 
bound is less than one, then, two-time- scale behavior is implied for any energy climb that the aircraft is 
allowed to perform. If however the resulting value of the upper bound turns out to be greater than one, 
then no conclusion can be drawn. The way to proceed in this latter case would be to evaluate a less 
conservative upper bound for e and apply the same reasoning. As it turns out, the less conservative the 
upper bound, the more work one has to perform in order to evaluate it. If all the upper bounds for e, 
evaluated for the entire envelope fail to yield any conclusions, the reasonable thing to do next is to 
evaluate e as a function of the energy E using all the assumptions made in the evaluation of reduced 
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solutions in aircraft energy climbs ('j fc 0,L=mg etc.). By evaluating in this sense, and at each energy level 
the absolute maximum value of the longitudinal loading factor we obtain a curve C on the e-E plane. The 
interesting properties of this curve are that for a given aircraft it need only be constructed once and that it 
lies above all other curves that may be evaluated similarly, but along the reduced solutions corresponding 
to specific problems. In other words, points on curve C represent upper bounds for e at the 
corresponding energy levels. The portions therefore of curve C where £ is less than one immediately 
show the energy levels where two-time-scale behavior (boundary layer transitions along constant E) can 
be expected. If there are any portions of curve C where e is greater than one, then no conclusions can be 
drawn as to the possible two-time-scale behavior at the corresponding energy levels. In the latter case 
one has again to evaluate a less conservative upper bound for e at these energy levels. Such less and less 
conservative upper bounds would of course eventually lead to the maximum value of the longitudinal 
loading factor evaluated as a function of E along the reduced solution corresponding to a specific 
problem. 

It should be clear now that if we are interested in the possible two-time- scale behavior of a vehicle 
along a particular trajectory (corresponding to a specific problem) then the least conservative upper 
bound for e would be the maximum longitudinal loading factor encountered along that (exact) trajectory. 
Calculating this upper bound would not be very useful since it would require the actual computation of 
the trajectory first. The idea presented above suggests that there may be a hope of avoiding this by 
starting with a more conservative upper bound, and proceeding with less and less conservative upper 
bounds. 

In order to demonstrate the above ideas in practice, numerical evaluations of e are presented in Figs. 
1-8 for four types of vehicles. For each type there is one plot showing the variation of the maximum 
longitudinal loading factor of the vehicle with energy E, and one or more plots showing the variation of 
the longitudinal loading factor with E along the reduced solution corresponding to a specific optimization 
problem. 

Figs. 1 and 2 show the results for an F-8 fighter 10 . The two optimization problems considered for 
this case were minimum time to a specified energy and minimum time to a specified downrange position. 
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The reduced solutions corresponding to these problems are obtained by maximizing (with respect to V) 
at each energy level the quantities (T-D)V for the former and [(T-D)V]/(Vo-V) for the latter, where V 0 is 
the maximum possible cruising speed of the aircraft and D is calculated at L=mg. Fig. 1 shows the actual 
paths in the envelope corresponding to these reduced solutions and to the maximum longitudinal loading 
factor of F-8. Fig. 2 shows the results for e evaluated along these climb paths. Since the maximum 
longitudinal loading factor of F-8 stays always below one in Fig. 2, it is reasonable to assume that for a 
any optimization problem, if the required energy gain is sufficient, the transitions to the reduced solution 
' will take place at nearly constant E, exhibiting the well known boundary layer structure. 

Figs. 3 and 4 show similar results for an F-15 fighter 11 . The two optimization problems considered 
in this case were again minimum time to a specified energy and minimum time to a specified downrange 
position. A maximum dynamic pressure constraint of 1500 lbf per square feet is imposed on the climb 
paths for this case. Due to the large thrust to weight ratio of F-15, the e levels in Fig. 4 are much higher 
than the ones corresponding to F-8 (compare with Fig. 2). In particular there is a small region at low 
energy where e exceeds one, implying that two time scale separation at these energy levels may not be 
appropriate for the above two optimization problems. 

Figs. 5 and 6 show the results for a conventional transport 12 . In this case however, the two 
optimization problems considered were minimum fuel to a specified energy and minimum fuel to a 
specified downrange position. The reduced solutions corresponding to these problems are obtained by 
maximizing (with respect to V and rj) at each energy level the quantities [(T-D)V]/f for the former and 
[(T-D)V]/(fVo-f 0 V) for the latter, where V 0 is the most fuel efficient cruising speed of the aircraft and f 0 
is the fuel consumption rate at this cruising flight condition 7 . It is interesting to note the very low levels 
of e in Fig. 6 (compare with Figs. 2 and 4), suggesting two time scale behavior for the entire envelope. 

Finally, Figs. 7 and 8 show the results for a hypersonic vehicle model, used by NASA and called 
“the Langley Accelerator” 13 . The only optimization problem considered in this case was minimum fuel to 
a specified energy, the reduced solution corresponding to which is obtained by maximizing the quantity 
[(T)T-D)V]/(mf) at each energy level (mass is not constant in this case). A maximum dynamic pressure 
constraint of 2000 lbf per square feet is imposed on the climb paths for this case. This particular vehicle 
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model employs a multimode propulsion system, sized for acceleration to orbital velocity and consisting 
of turbojet, ramjet, scramjet, and rocket cycles. Note that Figs. 7 and 8 correspond to a nonoptimal 
switching between the different propulsion cycles. Specifically, we first assume allowable operating 
ranges for the different propulsion cycles (expressed as bounds on the Mach number), and then, at each 
energy level we pick the cycle that maximizes [(ijT-D)V]/(mf) and (T\T-D)/m. The switching therefore 
from one cycle of propulsion to the other is abrupt, with no overlap. The actual points of cycle 
transitions are also shown in the figures. Note also that in Fig. 8 e is plotted against the speeds at which 

the constant energy contours intersect the zero altitude axis. The reason for this is that E is negative in 
this case. The actual value of e is likely to be much lower if a practical method for cycle transition is 
employed. Note finally that as the energy levels get higher and higher we approach the boundary of the 
envelope and e goes to zero. This is basically a characteristic of all aircraft (see also Figs. 2, 4, and 6) 
suggesting that the transitions to the reduced solution can be treated as boundary layers more succesfully 
at high energy levels than at low ones. The physical explanation for this comes from the behavior of the 
difference between the thrust T and the drag D. At low energy levels both the speed and altitude are low, 
implying that the thrust is high and the drag is low, so that T-D is high and can be used to change the 
energy during a transition. At high energy levels on the other hand either the speed or the altitude or both 
are high, implying that the difference T-D is low. Thus, transitions to the reduced solution at high energy 
levels can be expected to occur more or less by interchanging kinetic for potential energy (or vice versa), 
with the total energy staying nearly constant. 


Conclusions 

For both the conventional (subsonic-supersonic, flat Earth) and the trartsatmospheric (hypersonic, 
spherical Earth) flight regimes a systematic procedure was introduced to identify naturally a singular 
perturbation parameter e in the differential equations of motion. The procedure consists of using a set of 
arbitrary scaling constants to nondimensionalize all the variables of interest and then applying a set of 
conditions on these constants to put the resulting nondimensional equations of motion in the traditional 
singular perturbation format. Because the number of conditions is less than the number of constants, the 
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scaling constants cannot be uniquely specified. Thus, the resulting expression for e is in general quite 
arbitrary. There is, however, a particular choice of the scaling constants for which the values of the 
resulting e can serve as a hint for a possible two-time-scale behavior of the aircraft in question. The 
primary result of the paper is the demonstration that for this "useful" choice of the scaling constants the 
resulting £ is always equal to the maximum longitudinal loading factor of the vehicle, measured in units 
of g's. Two time scale behavior is suggested according to whether £ is less than one or not. Based on 
this result it is straightforward to see why singular perturbation methods applied to aircraft performance 
optimization problems have worked so well in the past. The maximum longitudinal loading factors 
associated with the majority of conventional aircraft are either less than one because the aircraft lack very 
high thrusting capabilities, or because they are restricted to be so for other reasons (structural, 
comforting, etc.). A few notable exceptions do occur for some modem fighters. This directly suggests 
that most conventional aircraft can be expected to exhibit two-time-scale behavior for almost any energy 
climb that they are allowed to perform. This then appears to be the reason for the past success of so 
many singular perturbation treatments of aircraft energy climbs. The implication for transatmospheric 
vehicles is rather direct. If we consider such a vehicle as a passenger transport, then, in order to assure 
passenger comfort it is only natural to impose as a constraint a maximum longitudinal loading factor for 
the vehicle that is less than one. Our work suggests that such a constraint would imply two-time-scale 
behavior for any type of energy climb that such a vehicle would be allowed to perform. Therefore, 
singular perturbation formulations of such maneuvers still appear to be promising, even with all the 
added complexities that the flight regimes of such vehicles can entail. 
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Fig. 1 Energy climb paths for an F-8 aircraft 
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Energy per unit weight (km) 
Fig. 2 Evaluation of 8(E) for an F-8 aircraft 
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Fig. 3 Energy climb paths for an F-15 aircraft 
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Fig. 4 Evaluation of e(E) for an F-15 aircraft 
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Fig. 5 Energy climb paths for a short-haul transport 
aircraft 
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epsilon 



Fig. 6 Evaluation of £(E) for a short-haul transport 
aircraft 
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Fig. 7 Energy climb paths for a generic hypersonic 
vehicle 
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Turbojet 



Fig. 8 Evaluation of 8(E) for a generic hypersonic 
vehicle 
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